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Using  the  method  of  embedding  a  given  equation  into  a  ^higher  dimensional 
problem4^"” examples  of  multiple  limit  point  solutions  emanating  from  a  single 
limit  point  are  constructed  within  the  context  of  algebraic  and  ordinary 
differential  equations.  In  both  instances,  the  linearized  problem  around  the 
limit  point  is  assumed  to  lead  to  non-self-adjoint  operators. 
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SIGNIFICANCE  AND  EXPLANATION 


In  order  to  solve  equations  of  the  form 

G(u,X)  -  0 

at  a  limit  point  where  u  =  u(X)  ceases  to  be  a  single  valued  function  of 
X,  and  du/dX  becomes  unbounded,  it  has  been  a  common  practice  in  recent 
years  to  imbed  the  equation  G(u,X)  =  0  into  a  'higher  dimensional1  problem 
so  that  the  Jacobian  of  the  enlarged  system  is  non-singular  at  the  limit 
point.  This  method  permits  one  to  seek  single  valued  solutions  for  u  =  u(e) 
and  X  =  X(e)  in  terms  of  a  new  parameter  e.  Despite  this,  it  is  common  to 
speak  of  the  solution  curve  passing  through  the  limit  point  as  exhibiting 
limit  point  bifurcation,  because  u  is  a  double-valued  function  of  X. 

Limit  point  bifurcations  arise  in  numerous  areas  of  mechanics,  such  as 
the  buckling  of  a  shallow  clamped  shell  under  uniform  pressure,  the  flow  of  a 
non-Newtonian  fluid  in  between  rotating  disks  and  in  chemical  reactor 
problems . 

In  the  present  paper,  the  possibility  of  two  curves  passing  through  a 
single  limit  point  is  discussed  and  is  shown  to  depend  on  the  dimension  of  the 
null  space  of  the  Jacobian  of  G(u,X)  at  the  limit  point,  as  well  as  the 
dimension  of  the  null  space  of  the  adjoint  of  the  Jacobian,  with  this  method 
it  becomes  possible  to  study  the  existence  of  multiple  limit  point  bifurca¬ 
tions  in  mechanics  and  engineering.  Finally,  because  the  theory  is  developed 
to  be  applicable  to  abstract  equations,  it  is  cast  in  terms  of  the  'Jacobian' 
being  a  Fredholm  operator  and  its  index  is  non-zero  because  the  dimensions  of 
the  null  spaces  of  the  operator  and  its  adjoint  are  not  equal. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
tunury  lies  with  MFC,  and  not  with  the  author  of  this  report. 


MULTIPLE  LIMIT  POINT  BIFURCATION  FOR 
FREDHPLM  OPERATORS  OF  NON-ZERO  INDEX 


R.  R.  Huilgol 

1 .  INTRODUCTION 

Let  X  end  Y  be  real  Banach  spaces  and  X  e  R  be  a  scalar.  Consider  a  nap 
G  :  X  *  R  ♦  Y  and  let  us  seek  solutions  of 

G(u,X)  -  0  (1) 

in  the  vicinity  of  a  known  singular  solution  (u°,X°)  e  X  x  R.  Assuming  G  to  be 

sufficiently  smooth,  we  shall  seek  a  solution  'arc"  (u (e),X(e)>,  with 

O  O  ,  |  o  o 

(u(0),X(0>)  3  (u  ,X  ) ,  depending  on  e  in  e  <  e  ,  in  the  neighborhood  of  (u  ,X  )# 

o 

If  we  expand  G(u,X)  about  (u°,X°)  we  find  that 

.  o  o  o  .  ,  o  1  o  o  o 

G(u,X)  -  G(u-u)+G(X-X)+  —  G  (u-u  Hu  -  u  )  + 
u  2  uu 

(2) 

+  G°^(u  -  u°)(X  -  X°)  +  -j  G°^(X  -  X°)(X  -  X°>  +  higher  order  terms 
o  3g  o.o 

In  ( 2 ) f  G  =  r—  (u  ,X  ),  and  the  other  operators  are  similarly  defined, 
u  rfu 

Let  G°  be  a  Fredholm  operator1  with  a  null  space  N (G°) ,  of  dimension  m  >  2,  and 

G?  be  such  that  it  does  not  lie  in  the  range  f?(G°)  of  G°,  i.e., 

X  u  u 

G°  $  R(G°)  .  (3.) 

X  T  u 

Now,  because  G°  is  Fredholm, 

X  -  N(G°)  •  X,  ,  (4) 

Y  -  R(G°)  •  Y,  ,  (5) 
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If  is  bounded,  the  argument  that  follows  is  unaffected.  If  is  unbounded,  its 

domain  which  is  dense  in  X  can  be  turned  into  a  Banach  space  via  the  graph  norm,  and  one 
can  proceed  with  requisite  modifications.  Details  are  left  to  the  reader. 
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where  Yy  has  the  same  dimension  es  that  of  the  null  space  of  the  adjoint  (G°)*  of  G°. 
If  we  seek  small  solutions  to  G(u,X)  -  0  when  G(u,X)  has  the  expansion  (2),  it  is 


clear  that  we  are  led  to: 


(u  -  u°)(e)  »  e*  +  e2w(c)  , 


(X  -  X°)(c)  -  e2c(e)  , 


where  ♦  e  W(G°),  w  e  ,  t  e  R.  As  (7)  implies,  on  every  solution  arc  bifurcating  from 

(u°,X°),  dX/de  -  0  at  E  -  0  [1,2]. 

Following  Keller  [1],  and  Decker  and  Keller  [2],  Hull go 1  [3]  showed  that  the  solutions 

o 

(6)  -  (7)  are  of  the  "multiple  limit  point  bifurcation"  type  provided  G^  is  bounded  and 

R(G°>  C  R(G°)  -  Y  ,  (8) 


R  <GX )  -  Yt  .  ( 

o 

Now,  a  glance  at  equations  (2)  -  (3)  will  show  that  because  G^  ■  R  ♦  y  ,  one  has 
dim  Y.|  •  1  (  moreover,  dim  Y^  »  dim  W((G°)*)  and  thus  the  Fredholm  index  of  G°  is 
m  -  1  y  0,  since  m  >  2. 


In  order  to  establish  (€)  -  (7),  one  can  fix  4  e  N(Gu>  and  define  a  new  map 


g  l  (  »  *  1  ♦  Y  through 


g(E,v,5)  - 


^  G(u°  +  e4  +  ev,  X°  +  e(),  e  y  o  , 


G°(4  +  v)  +  g°5,  e  -  o 


Clearly,  g(0,0,0)  »  0  because  ♦  e  NIG^) . 

Since  (8)  holds,  by  the  implicit  function  theorem  [4],  the  equation  g(e,v,5)  *  0  has 
a  unique  solution  v  -  v(E),  5  »  5(E),  for  |e |  <  Eq,  such  that  v(0)  -  0,  5(0)  ”  0. 

Using  the  smoothness  of  G,  one  can  write  v(c)  -  ew(e),  5(e)  -  e((e)  and  recover  (6)  - 


Now,  there  are  a  number  of  situations  as  in  ttopf  bifurcation  (Weber  [5]),  or  in  the 
example  below  when  it  is  computationally  attractive  to  add  an  extra  condition  to  (10)  - 
(11).  For  instance,  if  X  ■  rP,  we  know  that  we  can  choose  v  to  be  orthogonal  to  4* 


-2- 


More  generally,  if  G°  is  Fredholm  then  there  is  always  a  linear  operator  (e.g.,  a 
projection  operator)  C *  :  X  ♦  rF1  such  that  C  :  X1  +  (0)  and  C*  is  an  isomorphism 
between  N(G°)  and  rF1,  since  dim  N<G°)  ••  m. 

To  incorporate  such  a  C*,  it  is  essential  to  seek  the  solution  of  G(u,X)  »  0  about 
(u°,X°)  by  imbedding  the  original  problem  into  a  "higher  dimensional"  one.  That  is  to 
say,  we  solve  a  second  equation  q^(e,v,F.)  ”  0,  where  g1  :  R  *  X1  *  R  ♦  Y,  and  g,  is 
defined,  for  a  fixed  ♦  e  N(G°) ,  through 


91  (  e,v,C )  -  < 


1  G(u°  +  t*  +  ev,  X°  +  et),  £  0 

c*v,  £  f  0  , 


|G°<*  +  v)  +  G°5,  e 
|c*V,  E  -  0  . 


(12) 

(13) 

(14) 

(15) 


Clearly,  g1 (0,0,0)  «  0  and  in  order  to  employ  the  implicit  function  theorem  so  that 
9^(6, v,C)  ■  0  has  a  unique  solution  (v(e),£(e))  with  v(0)  •  0,  £(0)  »  0,  one  needs  to 
establish  that  the  "linear  operator" 


L' 


As 


(16) 


is  non-singular#  Note  that  A  arises  because  one  can  write  (14)  -  (15)  symbolically  as: 

v 

gl  (0,v,£)  -A  •  (17) 

It  was  proved  in  Theorem  I  of  (3]  that  A  is  non-singular,  provided 
(i)  R(G°)  ®  R(G°)  -  Y  , 

(ii)  N(C* )  -  Xx,R(c*)  -  R*  , 

(ill)  dim  N(G°)  -  m  , 

(iv)  W(C*)  n  N(G°)  -  {0}  . 

The  above  idea  has  obvious  extensions  to  the  case  where  X  e  in  (1)  -  see  nieorem  II  of 
[3].  This  is  not  pursued  here. 

Now,  on  putting  v(e)  “  Ew(e),  F.(e)  -  ec(e)  in  (6)  -  (7),  we  can  write  G(u,X)  *  0 
in  (2)  as: 


-3- 


r 


a 

G  G. 

- 

1  G°  *♦  +  0(e) 

w 

U  A 

w 

2  uu 

s 

* 

m  • 

C 

C  0 

C 

0 

-  - 

.  j 

- 

„  - 

(18) 


where  (6)  and  (7)  are  usad.  Since  dim  N(G°)  >  2,  the  possibility  exists  that  more  than 
one  solution  pair  (w,C)  can  be  found  by  solving  (18)  when  we  let  ♦  be  a  different 
linear  combination  of  {+\. . .  ,  $*)  which  span  N(G°). 

In  this  paper,  using  this  idea,  we  construct  two  examples  of  multiple  limit  point 
bifurcation,  when  X  e  R. 

2.  ALGEFRAIC  EQUATIONS 
Firstly,  let  us  consider  the  equations 

3  n  i  3 

Gi(S'X)  “  l  "  “41  *  *,<*  -  X°)  *  -  l  B  (u  -  u°)(u  -  u°) 

J-1  3  3  3  j,V-1  3  3  3  K  k  ( 19) 


+  higher  order  terms  -  0,  1  -  1,2  , 


where  is  the  2x3  matrix: 


11 


3G. 

_ i  ,  o  .o. 

(a  *x  >  ' 


(20) 


X  e  R.  *1, 

i  -  1,2 

,  are 

the 

conponents 

of  the 

vector  J  spanning  N(a  ) 

with  a  being  | 

the  adjoint 

of 

a. 

Clearly, 

4  t  R<s> 

.  Now,  let 

a  have  a  null  space  of 

dimension  2,  j 

.1 

.2 

* 

1 

spanned  by 

t 

and 

t 

respectively. 

Then, 

di»  N(a  )  -  1.  Now  we 

can 

rewrite  (19),  by  j 

appealing  to  (18), 

as: 

~ 

r  ** 

3  ^ 

] 

°11 

°12 

°13 

I 

wi 

1 

2 

l  Biik*A  +  0U) 

J,k-1  3*  3  * 

j 

°21 

°22 

°23 

1  *2 
I 

w2 

— — 

| 

-  - 

1 

2 

3 

l  <WA  +  tJ(e) 

j,k-1  *3*  3  * 

(21) 

♦; 

I  o 

1 

W3 

0 

m 

*1 

♦32 

1 

1  0 

1  J 

C 

w  - 

0 

-4- 


In  (21),  on  the  left  aide,  we  have  employed  the  fact  that  jfc1  and  are  orthogonal  to 

w.  Thus  the  lower  2*3  matrix  represents  C*.  On  the  right  side  of  (21),  we  have  taken 

1  2 

£  e  N(a)  to  stand  for  £  or  £  or  a  linear  combination  of  both. 

Since  rank  a  -  1,  the  row  [o^  a12  a13]  is  a  multiple  of  [a21  ct22  <*2;J),  with  a 
factor  Y,  say.  Provided  1(8^  “  Y82  jk  ^k  *  **  0  we  g8t  *  solution  C(0)  0.  This 


follows  from  the  fact  that  rank 


*11  *12  *13  *H 
°21  °22  °23  *2| 


■  2  and  thus  f  Y<»2. 


Now  we  can  have  two  distinct  solutions  for  C(0),  leading  to  multiple  limit  point 


bifurcation,  provided  (say) 
3 


As  ii  specific  example  let 


Then  N(a)  is  spanned  by: 


ij 


1  1 
1  1 


1 

-1 
0 

and  the  domain  X  of  a  has  the  decomposition: 


(22) 


(23) 


(24) 


X  =  K3  -  N(a>  •  X, 


X1  ”  span 


(25) 


3  2 

Similarly,  since  a  :  R  ♦  R  ,  the  range  space  Y  has  the  decomposition. 


s  R2  -  R(a 


)  •  Y 


where  R(a)  and  Y1  are  described  through 

f:l)  •  -  ■ 


i  ' 

i 

- 1 


,  dim  Y.  - 


(26) 


R(g)  “  span 


span 


1 


(27) 


Moreover,  £ 


[3 


■pens  N(a  )  as  well  and  does  not  lie  in  R(a ) . 


For  simplicity,  let  us  choose 


?111  "  1#  ®213  "  2' 


other  "  0 


(26) 


Then ,  if  we  take  £  ”  £  in  ( 2 1 )  we  get  at  e  »  0  s 


- 

r 

■  - 

1111 

i 

!  1 

2 

111-1 

k 

1 

0 

i 

i 

! 

■ 

1-10  0 

k 

0 

11-20 

ic° 

0 

_ 

v.  - 

_ _ 

(29) 


Now,  the  4*4  matrix  on  the  left  side  is  non-singular1  and  we  get  the  unique 
solutions : 

„°  -  w°  .  w°  -  -  1 

12  3  12 


(30) 


1 


Thus,  in  this  case,  the  solution  to  (19)  is  of  the  form 

(u  -  u°)(€)  ■  EjJ1  +  e2w°  +  0(e3),> 


(X  -  \°)(e)  -  eV  +  0(e3K 


(31) 


where  w°,  i  *  1,2,3,  and  C°  are  given  by  (30).  Had  we  employed  £2  in  (21),  we  would 

1  2 

get  a  solution  set  similar  to  (31),  except  that  £  is  replaced  by  £  and  the  new  values 


The  usefulness  of  incorporating  C*  in  (12)  -  (15)  becomes  clear  here.  For,  without  it, 
we  have  2  equations  in  (19)  for  Wj,  w2,  w3  and  (.  Of  course,  the  first  three  are  inter¬ 


related.  The  operator  C  makes  this  explicit  as  well  as  leading  to  a  square  matrix  A 
which  is  invertible. 


of  w°  and  c°  are  i 

(32) 

Hence  we  have  constructed  a  pair  of  algebraic  equations  yielding  oultiple  limit  point 
bifurcating  solutions. 

Now,  Decker  and  Keller  [2],  while  assuming  X  e  R,  require  to  have  a  Fredholm 

index  zero.  Therefore,  their  theory  is  inapplicable  to  the  example  chosen  here  since 
dim  Ma)  f*  dim  W(a  ),  while  the  example  follows  from  the  theorem  proved  in  (3). 

It  seems  to  be  the  simplest  example  of  multiple  limit  point  bifurcation,  since  this 
type  of  bifurcation  requires  (2]  dimN(G°)  >  2. 


3.  A  DIFFERENTIAL  EQUATION 

Once  an  example  such  as  in  (19),  (23)  -  (28)  has  been  constructed,  it  is  clear  how  to 
generate  others.  To  be  specific,  let  u°  ”  u°(x)  belong  to  X  =  L2(0,2x).  Now,  consider 
the  differential  equation  for  u  6  Hi 

2 

G(u,X)  ■  — -  (u  -  u°)  +  (u  -  u°>  +  (X  -  1 )♦  + 

d*  (33) 

♦  f(u  -  u°,  \  -  1,  (u  -  u°))  -  0  . 


In  (33), 

desires. 

Clearly, 

Now, 


f( 0,0,0)  *  0  and  f  is  not  linear  in  its  arguments  and  is  as  smooth  as  one 
The  boundary  conditions  associated  with  (33)  are  taken  to  bei 

u(0)  -  u(2ir)  .  (34) 

(u°,1)  la  a  solution  of  (33)  -  (34)  if  u°  e  X  meets  (34). 
the  adjoint  boundary  conditions  arei 

u(0)  -  u(2x )  -  0,  u’(0)  -  u'(2w)  .  (35) 


Hence,  if  we  let 


with  the  boundary  conditions  as  given  by  (34),  then 


N(G°)  -  span  {sin  x,  cos  x)  . 


(36) 


-7- 


However,  the  adjoint  of  G°  has  the  form  (G°)* 
conditions,  and  hence 


with  (35)  as  its  boundary 


N((G°)*)  «  span  {sin  x}  (37) 

leading  to  the  index  of  G°  being  1.  Now,  while  the  domain  V(G )  of  G  contains 
u  e  X  meeting  (34)  we  have  yet  to  specify  the  range  R(G)  of  G.  To  this  end,  set 
Y  ■  X\fcos  x}>  then  we  put  R(G)  -  Y.  Finally,  let  the  function  <>  be: 

♦(x)  »  sin  x  ,  (38) 

It  is  now  clear  that  (33)  is  similar  to  (19)  and,  therefore,  multiple  limit  point  bifurca¬ 
tions  arise  from  (u°,1)  provided  (cf,  (22)): 

f<VX  -  1,  V  *  f(*2'X  '  ^  *2)  (39> 

where  +.,(x)  •  sin  x,  ^(x)  »  cos  x. 


4.  CONCLUDING  REMARKS 

As  mentioned  in  [2],  there  are  numerous  examples  of  limit  point  or  turr.  g  nt 
bifurcation  in  the  literature,  especially  in  mechanics.  In  this  literature,  it  is  said 
that  there  are  two  solutions  near  a  turning  point  because  in  the  usual  (lul,X)  diagram 
near  (u°,X°),  there  are  two  values  for  lul,  for  each  X  ?  1°.  In  the  approach  adopted 
here,  following  (2,31,  such  a  solution  is  assumed  to  represent  a  single  limit  point 
solution  for,  in  terms  of  e,  the  solution  curve  through  (u°,X°)  has  a  unique  repre¬ 
sentation.  TO  he  specific,  if  c  «  0  corresponds  to  (u°,X°),  then  the  values  of  c  <  0 
and  e  >  0  correspond  to  distinct  values  of  (u,X),  respectively,  around  (u°,X°)>  or, 
we  have  a  smooth  solution  arc  through  the  limit  point. 

What  has  been  done  in  the  present  paper  is  to  construct  examples  of  double  solution 
arcs  through  the  limit  point.  In  the  spirit  of  (2,3),  these  are  called  multiple  limit 
point  bifurcations,  and  are  believed  to  be  the  first  known  examples  of  this  kind  in  the 


literature 
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